Tunneling Dynamics of Bose-Einstein Condensates 
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We study tunneling dynamics of atomic pairs in Bose-Einstein condensates with Feshbach res- 
onances. It is shown that the tunneling of the atomic pairs depends on not only the tunneling 
coupling between the atomic condensate and the molecular condensate, but also the inter-atomic 
nonlinear interactions and the initial number of atoms in these condensates. It is found that in addi- 
tion to oscillating tunneling current between the atomic condensate and the molecular condensate, 
the nonlinear atomic-pair tunneling dynamics sustains a self-locked population imbalance: macro- 
scopic quantum self-trapping effect. Influence of decoherence induced by non-condensate atoms on 
tunneling dynamics is investigated. It is shown that decoherence suppresses atomic-pair tunneling. 
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The atomic Bose-Einstein condensates JljH offer new 
opportunities for studying quantum-degenerate fluids. 
All the essential properties of atomic Bose-Einstein con- 
densed systems are determined by the strength of the 
atomic interactions. In contrast with the situation of the 
traditional superfluids, the strength of the inter-particle 
interactions in the atomic Bose-Einstein condensate can 
vary over a wide range of values through changing exter- 
nal fields. Hence one can manipulate and control conden- 
sate properties by varying the strength of interactions. 
The Feshbach resonance approach Q is considered as 
an effective one to alter the inter-atomic interactions in 
Bose-Einstein condensates. The magnetic-field-induced 
Feshbach resonances in an atomic Bose condensate have 
already been observed experimentally Theoretical 
studies of the ultracold atoms with Feshbach resonances 
HQ showed that the interactions responsible for the Fes- 
hbach resonances produce a second condensate compo- 
nent, a molecular condensate, and predict that tunneling 
of atomic pairs occurs between the atomic condensate 
and the molecular condensate. Recently, the molecular 
Bose-Einstein condensate has been produced experimen- 
tally [g. 

The purpose of this paper is to study tunneling dy- 
namics of atomic pairs between the atomic condensate 
and the molecular condensate. We show that in addition 
to oscillating tunneling current between the atomic con- 
densate and the molecular condensate, the nonlinearity of 
the tunneling dynamics sustains a self-maintained pop- 
ulation imbalance: macroscopic quantum self-trapping 
effect (MQST). We also discuss the influence of decoher- 
ence induced by non-condensate atoms on tunneling dy- 
namics and find that decoherence suppresses atomic-pair 
tunneling. 

The binary atom Feshbach resonances are hyperfine- 
induced spin-flip processes that bring the colliding atoms 
to a bound molecular state of different spins, and then re- 
turn an unbound state. These processes can be described 
by Hamiltonian (% = 1) 



H 



FR 



a / drip+(r)ip a (r)ip a (r) +h.c, 



(1) 



where ip m (r), V>^(r) (^a r )> - t ^( r )) are the annihilation 
and creation field operators of the molecules (atoms), a 
stands for the coupling constant. The Hamiltonian Hfr 
together with atomic, molecular, and atom-molecule in- 
teraction Hamiltonians (K = 1) 



H a = 



H„ 



^ Q + (r)[-^:V 2 + y(r)]^(r) 

dri+(r)[- i i ? V 2 + V(v) + e]^ m (r) 
-^f [ dr^+(r)^+(r)V; m (r)^ m (r), 



A' / dr^+(r)^+(r)^ a (r)^ m (r), 



H„ 



forms a total Hamiltonian 



H = H n 



H 



FR, 



(2) 

(3) 
(4) 

(5) 



which governs the dynamics of the system under our con- 
sideration. Here V(r) represents the trapped potential, 
^'a(m) ~ ^ 7Ta a(m) /((2)Af) with M being the atomic mass 
and a a (m) the scattering length, A' denotes the coupling 
constant of atom-molecule interaction. The detuning e 
linearly depends on the magnetic field e oc B — Bq with 
Bq being the resonant magnetic field. 

For small atomic and molecular condensates the 
atomic and molecular field operators can be approxi- 
mated as V> Q (r) = acj) a (r), ip m (r) = b(f> b (r) where <p a (r) 
and </>b(r) are real normalized mode functions for the 
two condensates, and a and b are associated mode an- 
nihilation operators which satisfy the standard bosonic 
communtation relations. Then the total Hamiltonian be- 
comes the two-mode Hamiltonian 
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H = H a + H 1 , (6) 
H = u) a a + a + u) b b + b + X a a +2 a 2 + X b b +2 b 2 + Xa + ab + b, 
H' = a{b + a 2 + ba+ 2 ), 

where ab + a 2 describes the annihilation of a atomic pair 
in the atomic condensate and the creation of one molecule 
in the molecular condensate thereby transferring a pair of 
atoms from the atomic condensate to the molecular con- 
densate with a being the corresponding tunneling cou- 
pling constant. The Hermitian conjugate part aba +2 de- 
scribes the reverse process. Hence, what the Feshbach 
Hamiltonian H' describes is not only a three-body re- 
combination process of the molecular formation from a 
pair of atoms but also a tunneling process between the 
atomic condensate and the molecular condensate. 

In general, the Hamiltonian (6) can not be exactly 
solved, but it can be perturbatively solved in the off- 
resonant regime with a large detuning e 3> 0, by which we 
mean here that e greatly exceeds the Feshbach-resonant 
interaction energy, so we can treat H' as a perturbation. 

Let Ho\n,m) = En,ln\n, m) where \n,m) is an eigen- 
state of the number operators a + a and b + b defined by 
a + a\n, m) = n\n, m) and b + b\n, m) = m\n, m). It is easy 
to find that 



E n}m — nUJ a + mUJ b + (n 2 - n)X a 

+(m 2 — m)X b + nmX. 



(T) 



For simplicity, we consider the nondegenerate case and 
assume that H\ip n , rn ) = E n . m \ip niTn }. Then the pertur- 
bative energy and eigenstate are found to be 



E nL + a [ a n,mVn(n - l)(m + 1) 



+&n,ro\/ ( n + !)(«■ + 2)m] 



(8) 



\ipn,m) ~ A n , m [\n,m) + a„, m |n - 2,m+ 1) 

+b n>m \n + 2,m-l)}, (9) 

where the normalization constant A n>m and two coeffi- 
cients a n .m and b niTn are given by 



A — 



1 





+ b 2 


oi\J n(n - 


- l)(m+ 1) 


p(0) 


P (0) 


U\] (n + 


l)(n + 2)m 



™' m ~ p(0) _ p(0) 



(10) 



Let the two condensates be initially in a state |^(0)) = 
En m Cn,m(0)\ipn.m), the time evolution of the wave 
function is then given by the expression 



|tt(t)) « > C n>m (Q)e- 1 ^ '|Vn, m ) 



n.m 



(ii) 



In order to investigate tunneling dynamics, we intro- 
duce the population difference 



P(t) = n«(t) - 2rc 6 (i), 



(12) 



where n a nA (t) is the number of atoms (molecules) in the 
atomic (molecular) condensate at time t. 

Now let us assume that initially both the atomic and 
molecular condensates are in a Fock state |m, 112). From 
Eqs. (11) and (12) we find that 
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P(t) = P (n 1 ,n 2 ) +^Pi(ni,n 2 )cos[w i (ni,n 2 )£], (13) 

i=l 

where the oscillation frequencies are given by 
wi(ni, n 2 ) = E niin2 — E ni +2,ri2-ii 

^2(^1, n 2 ) = E nii n 2 — i?„ 1 _2 1 n 2 + lj 

W3(ni,n 2 ) = £ ; „ 1+ 2,n 2 -l — Eni -2,n2+l J (14) 

and the coefficients Pj(ni,n 2 ) are defined by 

Po(ni,n 2 ) = A 2 ni n2 [{ ni - 2n 2 ) + {n x - 2n 2 - 4)c£ ljna 
+(n 1 -2n 2 +A)bl un2 ] 

+^n 1 +2,n 2 -l a n 1 +2,n 2 -l[(«l - 2n 2 + 4) 

+(ni - 2n 2 )a 2 ni+2 n2 _ 1 
+(n 1 -2n 2 + 8)b 2 ni+2jn2 _ 1 ] 

+^i-a,n a +i^i-a, na +i[(»*i - 2 "2 - 4) 
+(ni - 2n 2 - 6)a ni _ 2 n2+1 
+(n 1 -2n 2 )6, 2 li _ 2 ^ 2+1 ], (15) 

Pi(ni,n 2 ) 

x [(ni - 2n 2 )a„ 1+2 ,„ 2 _i 
+(ni - 2n 2 + 4)6 nii „ 2 ], 

P 2 (n 1 ,n 2 ) A ni — 2,712 + 1 ^Tll— 2,712 + 1 

x[(rai - 2n 2 - 4)a ni ,„ 2 
+(ni - 2n 2 )6 ni _ 2 ,„ 2+ i], 
PaC^i, ^2) = 2A„ 1+ 2 ! „ 2 -iA ni _ 2 ^ l2+ ia 2 i+2 „ 2 _ 1 

x^ 1 - 2 ,„ 2 +iK-2n 2 ). (16) 

From Eq. (12) we see that the population difference 
between the atomic condensate and the molecular con- 
densate exhibits oscillating behaviors. Especially, we can 
obtain the nonzero time average of the population differ- 
ence labeled by P given by 



P = Po(»i,n2), 



(17) 



which implies that there is a self-locked population imbal- 
ance between the atomic condensate and the molecular 
condensate. This is the MQST which occurs in the usual 
two Bose condensate system |l0|,[ll} as well. It is easy to 
check that the MQST vanishes when the nonlinearities in 
interactions are absent. Hence the MQST phenomenon 
is a nonlinear effect. 
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Eq. (12) indicates that the population imbalance be- 
tween the atomic condensate and the molecular conden- 
sate exhibits nonlinear oscillations with the time evolu- 
tion. It is these oscillations that leads to a Josephson-like 
tunneling current between the atomic condensate and 
molecular condensate, which can be defined as I(t) = 
P(t)/N with N being the total number of atoms in the 
atomic and molecular condensates. Making use of Eq. 
(13) it is easy to find that 



'(*) = -E 



P i {n 1 ,n 2 )uji{ni,n2) 
ni + 2n 2 



sin[wi(ni, n 2 )t]. (18) 



From Eqs. (7)-(16) and (18) we see that the tunnel- 
ing of the atomic pairs depends on not only the tun- 
neling coupling between the atomic condensate and the 
molecular condensate but also the inter-atomic nonlin- 
ear interactions and the initial number of atoms in these 
condensates. 

In order to further understand the influence of inter- 
atomic interactions on tunneling dymamics of the system 
under our consideration, let us specialize to the case of 
the atomic condensate initially being in a number state 
\N), and the molecular condensate initially being unpop- 
ulated, i.e., |^(0)) = \N, 0). In this case, the population 
difference and the tunneling current between the two con- 
densates given by the expressions 

P(t) = P (N,0) + P 2 (N, 0)cos[w 2 (N,0)t], (19) 
I{t) = -I am wa[u2(N,0)t], (20) 

where we have set I am = P 2 (N 7 0)uj 2 (N,0)/N. 

It is easy to see that the amplitude I am and the fre- 
quency co 2 (N, 0) depend upon the initial number of atoms 
in the condensates, the tunneling coupling a, the nonlin- 
ear interaction strengths A a , A&, and A. To find out how 
the initial number of atoms in the condensates and the 
interaction strengths affect the tunneling current, in Fig- 
ure 1 we plot the amplitude of the tunneling current as a 
function of interaction strengths for different initial num- 
ber of atoms when the initial state is |^(0)) = \N, 0), 
A a = Aj, = A, and uj a = 2u>i,. Figure 1 indicates that 
the amplitude of of the tunneling current is almost in- 
dependent of the initial number of atoms, the tunnel- 
ing coupling, and inter-atomic nonlinear interactions in 
the regime of weak (strong) tunneling coupling (nonlinear 
couplings) < a/X < 4. However, the tunneling coupling 
and the inter-atomic nonlinear interactions strongly af- 
fect the amplitude of the tunneling current in the regime 
of strong (weak) tunneling coupling (nonlinear couplings) 
a/X > 4. From Figure 1 we can see that the amplitude 
of the tunneling current increases with increasing both 
the initial number of atoms and the tunneling coupling 
in the regime of strong (weak) tunneling coupling (non- 
linear couplings) a/X > 4. 

In Figure 2, we display the scaled frequency of the 
tunneling current, lu/X~ lj 2 (N,0)/X, as a function of in- 
teraction strengths for different initial number of atoms 



when the initial state is |*(0)) = \N,0), X a = X b = A, 
and uj a = 2ujb- It is interesting to note that from 
Figure 2 we can see that there exists a zero-frequency 
point, labeled by D. From Eqs. (14) it is straight for- 
ward to see that the zero-frequency point is a degener- 
ate point of energy of the system under our consider- 
ation, at which the nondegenerate perturbation theory 
is broken. From Eqs. (7), (8), (14), and (20) we can 
find that the degenerate point is given by the expression 
a/X = [(3N-A)(7N- 2Q)/2(N - 2)(N - 3)] 1 / 2 . Figure 
2 indicates that on the left hand side of the degenerate 
point D the scaled frequency decreases with increasing 
the tunneling coupling and/or decreasing the nonlinear 
interaction strengths , and increases with increasing the 
initial number of atoms. On the other hand, on the left 
hand side of the degenerate point D the scaled frequency 
increases with increasing both the tunneling coupling and 
the initial number of atoms, and/or decreasing the non- 
linear interaction strengths . 

However, it is customary to consider a Bosc-Einstein 
condensate to be in a coherent state, associated with a 
macroscopic wave function with both amplitude and a 
phase, the presence of which is originated from Bose bro- 
ken symmetry. Assume that the two condensate are ini- 
tially in the coherent states \a) and \/3), which are eigen- 
statcs of a and b, repectively, then we have 



Cn,m(0) — 



-An ' 



. a n (3 r - 



cxp(H 2 + |/3| 2 ) L V^! 

y/(n- 2)!(m+ 1)! 
6„, m a" +2 /3"" 1 
+ v/(n + 2)!(m^T)r 



(21) 



Making use of Eqs. (11), (12) and (21) we find the 
expression of the population difference 

p (t) = 2 ^2{^Po(n, m)C(n, m; n, m) 



+Pi(n, m)C(n, m; n — 2, m + 1) 
x cos(.E„ iTO — £'„_2,m+l)i 
+p 2 (n, m)C(n, m; n + 2, m — 1) 
x cos{E n>m - E n+ 2 >m -i)t 
+p 3 (n, m)C(n, m; n + 4, m — 2) 
x cos(-E„ >m — £ , n+ 4 jTra _ 2 )i}, 



(22) 



where we have introduced the notations 

Po(n, m) — (n — 2m) + (n — 2m — 4)a n r 

+(n - 2m + 4)fe 2 jm , 
Pt(n, m) = (n-2m- 4)a„, m , 
p 2 (n, m) = (n - 2m + 4)6„, m , 
p 3 (n, m) = (n - 2m + 4)a n+4 , m -2 Vm, 



(23) 



(24) 
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and 

C{n,m-n',m r ) = C n , m (0)C* w (0). (25) 

From Eq. (22) we can get the nonzero time average of 
the population difference labeled by P given by 

P = y^C(n, m; n, m)po(n, m), (26) 

which implies that there exists the MQST between the 
atomic condensate and the molecular condensate. 

From Eq. (22) we can obtain the Josephson-like tun- 
neling current 

I(t) = -y^[ii(n, m) sin(_E„ jm - £ n _ 2 , m +i)t 

-7 2 (n,m) sm(E niTn - E n+2 , m -i)t 
-I 3 (n,m) sm(E niTn - E n+itm - 2 )i\, (27) 

with 

I\(n, m) — 2pi(n, m)C(n, m; n — 2, m+ 1) 
x(E 

n.rn 

I 2 {n, m) — 2p 2 (n, m)C(n, m; n + 2, m — 1) 
x(E 

n,m 

E n+2 , m -t)/N, (28) 
/3(n, m) — 2p^(n, m)C(n, m; n + 4, m — 2) 
x(E 

n.rn 

where 7V=|a| 2 + 2|/3| 2 is the total number of the atoms 
in the two condensates. 

We now discuss the effect of the decoherence. In exper- 
iments on trapped Bose condensates of atomic gases, con- 
densate atoms continuously interact with non-condensate 
atoms (environment). As is well known, interactions 
between a quantum system and environment cause two 
types of unwelcomed effects: dissipation and decoherence 
Jl2[ . The dissipation effect, which dissipates the energy 
of the quantum system into the environment, is charac- 
terized by the relaxation time scale r r . In contrast, the 
decoherence effect is much more insidious because the 
coherence information leaks out into the environment in 
another time scale r^, which is much shorter than r r . 
Since macroscopic quantum phenomena in Bosc-Einstein 
condensates mainly depend on rather than r r , the dis- 
cussions in present paper only focus on the decoherence 
problem rather than the dissipation effect. 

We use a reservoir consisting of an infinite set of har- 
monic oscillators to model environment of condensate 
atoms and molecules in a trap, and assume the total 
Hamiltonian |l3| to be 

Ht = H + Y. "Mk + F({S}) £ c k (b\ + b k ) 

k k 

+^({s}) 2 e4 ( 29 ) 

k W k 



where the second term is the Hamiltonian of the reser- 
voir. The last term is a renormalization term. The third 
term represents the interaction between the system and 
the reservoir with a coupling constant Cfc, where {S} is a 
set of linear operators of the system or their linear com- 
binations in the same picture as that of H, F({S}) is an 
operator function of {S}. In order to enable what the in- 
teraction between the system and environment describes 
is decoherence not dissipation, we require that the lin- 
ear operator S commutes with the the Hamiltonian of 
the system H. Then, the interaction term commutes 
with the Hamiltonian of the system. This implies that 
there is no energy transfer between the system and its 
environment. So that it does describe the decoherence. 
The concrete form of the function ^({5}), which may 
be considered as an experimentally determined quantity, 
may be different for different environment. 

The Hamiltonian Ht can be exactly solved by using 
the unitary transfor- 

mation U = exp[H ^ k {ck/oJk){b\ - bk)}- Corresponding 
to the Hamiltonian (29), the total density operator of 
the system plus reservoir can be expressed as py(t) — 

e -i5*^-i e -«E k u » 8 I 8 *^ r (0)C)'- 1 e < *^* w * S * 8 *&e i * t . 
We assume that the system and reservoir are initially in 
thermal equilibrium and uncorrelated, so that /3t(0) = 
/3(0) <g) pr, where p(0) is the initial density operator of 
the system, and pu the density operator of the reser- 
voir, which can be written as pr — Y\kf>k{^) with pk(0) 
is the density operator of the fc-th harmonic oscillator 
in thermal equilibrium. After taking the trace over the 
reservoir, we can get the reduced density operator of the 
system, denoted by p(t) — tr^pxit), whose matrix ele- 
ments in the eigenstate representation of H are explicitly 
written as 

P(rn',n')( m ,n)(t) = \P(m' ,n')(m,n) (0) |e~ 7 < m '. "')<-">") (t) 

Xe -'>( m ',n')(m.»)(t), (30) 

where the damping factor and the phase shift are defined 
by 

7(m' ,n'){m,n)[t) = v 2 _(m , n '; m, n)Q 2 (t) , (31) 
<i>( m >,n>)( m ,n){t) = v+(m' ,n' ;m,n)«_(m' ,n' ;m,n)Q 1 (t) 

(32) 

where we have introduced the following notations: 
v±(n,m;n',m') = F({S(n,m)}) ± F{{S{n' ,m')}), (33) 

P(m',n')(m,n)(ty = \P(m' ,n')(m,n) (0)|e~ l9 ( m '.™'>( m> ™) , (34) 

and the two reservoir-dependent functions are given by 

Q 1 (t)=[ duJ(u)— ^sin(wt), (35) 
Jo ^ 

Q 3 (t)=2 r^JH^!M sin 2 (^)coth(^). (36) 
Jo w 2 2 
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Here we have taken the continuum limit of the reservoir 
modes: ^ fe — > J °° du>J(w), where J(u>) is the spectral 
density of the reservoir, c(u>) is the continuum expres- 
sion for Cfe, and j3 = 1/fcgT with fee and T being the 
Boltzmann constant and temperature, respectively. 

Eq. (30) indicates that the interaction between the 
system and its environment induces a phase shift and a 
decaying factor in the reduced density operator of the 
system. We now consider the population difference be- 
tween the atomic condensate and the molecular con- 
densate in the presence of the decoherence defined by 
P(t) = Trp(t)(n a - 2n b ). Making use of Eq. (30), Wc 
find that 

n,m)(n,m) (0)1 

COS @(n,m)(n,m) 

+pi(n, m)|p(„ ;m )(„_ 2 , TO +i)(0)| 

X COS 0(„, m )(„_2, ro +l) (t) e -T(».m)(»-=".m+l) (t) 

+P2(n,m)|p ( „ ;ro )(„ +2;ro _i)(0)| 

xcos0 ( „ !m)(ll+2!m _ 1) (t) e -^"-)("+ 2 -— DM 

+P3(n,m)|p(„ ;ro )(„ +4)ro _2)(0)| 

x cos</» ( „ ;ro)( „ +4;ro _2)(i)e-T(".-Kn + 4, ro - 2 )(t)}. (37 ) 

Then, the tunneling current is given by 

I{t) = -^{pi(n,m)|p (niTO )( n _2, m+ i)(0)| 

x [7(n,m)(n-2,m+l)(i) + 0(n,m)(n-2,m+l) (*) 
X Sin 0( n , m )(„_2,m+l) (t)]e-^.^-2, m+ r, (*) 
+p 2 (™, "l)|P(n,m)(n+2,m-l)(0)| 

x [7(n,m)(n+2,m-l)(£) + 0(n,m)(n+2,m-l) (*) 
X Sin ^( n , m )( n+ 2, m -l) (t)] e -T(».m)(»+='.m-i) (*) 
+P3(n,m)|p ( „ ;m )(„ +4;m _2)(0)| 

x [i(n,m)(n+A,m-2){t) + 0( n ,m)(n+4,m-2) (*) 

x sin^ (nim)(n+4im _2)(t)]e-T<"-K«+ 4 -- 2 )W}. (38) 

From Eqs. (37) and (38) we can immediately draw 
one important qualitative conclusion: since 7( n ,m)(n',m') 
is positive definite, the existence of the decoherence is 
always to tend to suppress the population difference and 
tunneling current between the atomic condensate and the 
molecular condensate. 

From Eqs. (31), (32), and (35)-(38) we see that 
all necessary information about the effects of the en- 
vironment on the population difference and the tun- 
neling current is contained in the spectral density of 
the reservoir. To proceed further let us now special- 
ize to the Ohmic case with the spectral distribution 
J(uj) = [i]Ll! / c 2 (uj)] exp (— ui/u c ), where lo c is the high 
frequency cut-off, r\ is a positive characteristic param- 
eter of the reservoir. With this choice, at low tempera- 
ture the functions Q\(t) and Q 2 (i) are given by the ex- 
pressions Qi(t) = ?7tan -1 (a> c t) and Q 2 (i) = ^{^lnfl + 



(ui c t) 2 ] + m[J^ sinh(^)]}. In particular, At zero temper- 
ature and in the meaningful domain of time uj c t ^> 1 , wc 
have Qi{t) = i]/(uj c t 2 ), and Q2(t) = rjln(uj c t). Then wc 
find 

n,m)(n,m) (0)1 

COS @(n,m)(n,m) 

n :: m 

+pi(n,m)\ 

P(n,m)(n— 2,m+l) (0)1 
X COS 0(„, m )( n _2,m+l) {t){uj c t)~^nm) 2 

+p 2 (n, m)|p(„ )m )( n+ 2, m _i)(0)| 

X COS ^( n , m )( n+ 2, m -l) (t)(uj c ty (v ^) 2 

+p 3 (n, m)|p (n!m)(n+ 4 !m _ 2 )(0)| 

x cos^ (nim)(n+4im _2)(t)(a; c t)-( 7 '3"n m ) 2 } ) (39) 

n.m 

+ ^ sin (nm)(n _ 2m+1) (t)] M )-K-„ m ) 2 

UJ c t 

+M n > m )\P(n,m)(n+2, m -i)(0)\[(V2nm.) 2 

+ ^sin0 ( „ !m)(n+2 , m _ 1) (t)]M)-( f '-™) 2 

+p 3 (n,m)|p ( „ !m)(rl+4!m _ 2) (0)|[(77 3nm ) 2 

+ sin0 ( „ !m)(n+4 , m _2)(t)]M)- (? '-"" )2 }, (40) 

U! c 

where we have used rif nm = y/rjv±(n,m;n — 2, to + I), 
vf„m = Vvv±(n,m;n+2,m-l),r)f nm = ^qv ± (n,m;n+ 

4, to- 2), and 77+" = ptumVinm- Eqs- ( 39 ) and ( 40 ) in- 
dicates that the tunneling current decays away according 
to the "power law" , where we have noted that the decay- 
ing factors can not be taken outside the summation on 
the r.h.s. of Eqs. (39) and (40). 

In summary, we have studied tunneling dynamics of 
atomic pairs in atomic Bose-Einstein condensates with 
Fcshbach resonances, and shown that the tunneling of 
the atomic pairs depends on not only the tunneling cou- 
pling between the atomic condensate and the molecu- 
lar condensate, but also the inter-atomic nonlinear inter- 
actions and the initial number of atoms in these con- 
densates, especially, we have shown that the tunnel- 
ing coupling and the inter-atomic nonlinear interactions 
strongly affect the tunneling of atomic pairs in the regime 
of strong (weak) tunneling coupling (nonlinear couplings) 
when the atomic condensate is in a number state and the 
molecular condensate in the vacuum state. This implies 
that the tunneling of atomic pairs between the atomic 
condensate and the molecular condensate can be manip- 
ulated and controlled by varying the tunneling coupling 
and/or inter-atomic nonlinear interaction strengths . We 
have revealed the existence of the MQST between the 
atomic condensate and the molecular condensate. The 
MQST is a kind of nonlinear effects which vanishes in the 
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absence of the inter- atomic nonlinear interactions. We 
have also discussed the influence of decoherence induced 
by non-condensate atoms on the tunneling dynamics, and 
shown that decoherence suppresses the atomic-pair tun- 
neling. Finally, it should be mentioned that inelastic col- 
lisions between the atomic and molecular conden- 
sates and reservoir may affect the atomic-pair tunneling 
between the atomic and molecular condensates. Influ- 
ence of inelastic collisions between the atomic and molec- 
ular condensates can be taken account into through in- 
troducing an imaginary part in the interaction strengths 
A and Af, |(|. Inelastic collisions between the system of 
the atomic-molecular condensates and reservoir lead to 
dissipation . A detailed investigation on the dissipation 
problem of the system of the atomic-molecular conden- 
sates is beyond the scope of the present paper, and will 
be given elsewhere. 
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FIG. 1. Amplitude of the tunneling current as a func- 
tion of interaction strengths for different initial num- 
ber of atoms when the initial state is 1^(0)) = \N, 0), 
A a = A5 = A, and uj a = 2cj . Here we have set 

LO = LJ 2 {N,0). 

FIG. 2. Scaled frequency of the tunneling current as a 
function of interaction strengths for different initial num- 
ber of atoms when the initial state is |^(0)) = \N, 0), 
A a = A;, = A, and uj a — 2cj . Here we have set 

LU = UJ 2 (N,0). 
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